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•  The  model  is  used  for  cell-level  design  analysis. 
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A  2D  direct  borohydride  fuel  cell  (DBFC)  model  has  been  developed  to  explore  the  prospective  perfor¬ 
mance  of  this  technology,  for  a  cell  with  fast  selective  electrocatalysts  and  a  selective  membrane.  In  the 
modeled  DBFC,  a  Nation  membrane  in  the  Na+  form  separates  flow  channels  with  aqueous  fuel  (0.1 
-0.5  M  NaBH4/4  M  NaOH)  and  oxidizer  (4  M  H2O2/4  M  H2SO4).  Electrochemical  reactions  occur  on 
catalyst-coated  channel  walls.  The  electrocatalysts  are  selective  for  complete  BH4  oxidation  and  H2O2 
reduction,  the  reactions  have  fast  forward  rate  constants,  and  only  Na+  and  H2O  cross  the  membrane.  The 
model  captures  interfacial  charge  transfer  reactions  and  complex  transport  in  the  flow  channels  and 
membrane.  Results  show  that  current  density  and  voltage  efficiency  vary  by  >50%  from  inlet  to  outlet 
due  to  concentration  boundary  layer  development.  The  BH4  concentration  boundary  layer  limits  peak 
power  density,  despite  migration  and  fuel  utilizations  below  10%.  Power  density  increases  with  BH4  inlet 
concentration  and  fuel  flow  rate,  but  at  the  expense  of  lower  fuel  utilization.  Water  crosses  the  mem¬ 
brane  up  to  14  times  its  production  rate  at  the  anode.  Low  fuel  utilization  and  water  imbalance  suggest 
the  importance  of  system  designs  with  reactant  recirculation  and  water  recovery. 

Published  by  Elsevier  B.V. 


1.  Introduction 

The  direct  borohydride  fuel  cell  (DBFC)  configuration  in  this 
study,  illustrated  in  Fig.  1,  uses  a  Nation  membrane  in  the  Na+  form 
to  separate  flow  channels  containing  aqueous  fuel  (0.1— 0.5  M 
NaBH4/4  M  NaOH)  and  oxidizer  (4  M  H2O2/4  M  H2SO4)  solutions. 
Borohydride  (BH4)  is  oxidized  on  the  anode  catalyst  at  the  fuel 
channel  wall  and  hydrogen  peroxide  (H2O2)  is  reduced  on  the 
cathode  catalyst  at  the  oxidizer  channel  wall.  The  catalyst-coated 
channel  walls  function  as  current  collectors.  DBFC  experiments 
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using  alkaline  NaBH4  and  acidic  H2O2  solutions  with  Nafion 
membranes  [1—3]  and  similar  separated  electrode  cell  topologies 
[1,4—6]  have  been  reported  in  the  literature. 

In  an  ideal  DBFC,  complete  BH4  oxidation  (Reaction  (R1 ))  yields 
8e~  with  a  standard  reduction  potential  of  f-Vi,  and  complete  H2O2 
reduction  in  acidic  media  (Reaction  (R2))  yields  2e~  with  a  standard 
reduction  potential  of  £°r2  [7j. 

B02  (aq)  +  6H20  +  8e“  ±4  BH4  (aq)  +  80H~ 

n  V  R 1  ) 

Eg,  =  -1.24  V  vs.  RHE 

H202(aq)  +  2H+  +  2e~  ;=±2H20  Eg2  =  1.763  V  vs.  RHE 

(R2) 

The  electrochemical  circuit  is  closed  by  sodium  (Na+)  crossing 
the  membrane,  which  also  drives  a  net  water  flux  from  fuel  solution 
to  oxidizer  solution  via  electro-osmotic  drag.  The  net  reaction 


R.O.  Stroman,  CS.  Jackson  /  Journal  of  Power  Sources  247  (2014)  756-  769 


757 


(Reaction  (R3))  consumes  BH4  and  CII-T  at  the  anode  and  H2O2  and 
H+  at  the  cathode  to  produce  B02  and  H20,  with  a  theoretical  cell 
potential  of  £°r3. 

BH4  (aq)  +  80H~  +  4H202  +  8H+-+B02  (aq)  +  14H20 
fg3  =  3.01  V 

(R3) 

This  DBFC  differs  from  common  low  temperature  fuel  cells  such 
as  the  proton  exchange  membrane  fuel  cell  (PEMFC)  and  direct 
methanol  fuel  cell  (DMFC).  The  catalyst  in  PEMFCs  and  DMFCs  is  on 
the  membrane  to  enable  H+  participation  in  the  electrochemical 
reactions,  but  the  catalyst  in  this  DBFC  can  be  elsewhere  because 
the  aqueous  electrolytes  in  each  channel  support  ion  transport. 
Channel  transport  in  PEMFCs  and  DMFCs  is  governed  by  convection 
and  diffusion,  but  in  a  DBFC  migration  also  contributes.  Finally, 
slower  diffusion  in  the  liquid  phase  favors  development  of  steeper 
concentration  gradients  in  a  DBFC. 

Interest  in  DBFCs  has  grown  due  to  several  desirable  features  of 
this  technology.  The  reactants  can  be  stored  at  sufficiently  high 
concentrations  to  give  theoretical  energy  densities  greater  than 
state-of-the-art  rechargeable  batteries  [81.  The  aqueous  reactant 
solutions  can  be  stored  at  ambient  temperature  and  pressure.  The 
theoretical  cell  voltage  is  3.01  V  when  the  oxidizer  is  H202  in  acidic 
media,  which  substantially  exceeds  cell  voltages  of  H2— 02  fuel  cells. 
In  addition,  the  reactant  chemical  energy  can  be  converted  directly 
to  electricity  without  any  chemical  preprocessing,  thereby  reducing 
power  system  complexity  and  improving  reliability.  Despite  these 
advantages,  challenges  remain  before  DBFCs  can  realize  their 
theoretical  promise.  Continued  progress  toward  meeting  those 
challenges  has  promoted  interest  in  DBFCs  for  applications  where 
energy  density  is  important,  particularly  for  air-independent  pro¬ 
pulsion,  remote  sensors  and  portable  electronics. 

For  DBFC  systems  to  be  practical,  higher  performance  electro¬ 
catalysts  and  membranes  must  be  developed  [9,10],  Experiments 
reported  in  the  literature  show  poor  catalytic  activity  and  the  loss  of 
reactants  to  unproductive  side  reactions.  BH4  crossover  has  been 
reported,  leading  to  mixed  potentials  and  decreased  cell  voltage. 
Membrane  conductivity  suffers  due  to  the  low  mobility  of  large 
ions  (such  as  Na+)  in  common  membrane  materials. 

Reactions  (R1 )  through  (R3)  assume  complete  oxidation  of  BH4 
at  the  anode  and  reduction  of  H202  at  the  cathode.  In  practice,  there 
are  intermediate  steps  and  competing  reactions  which  can  occur  at 
each  electrode.  At  the  anode,  the  complex  network  of  charge- 


transfer  steps  and  other  non-electrochemical  catalytic  reactions 
that  constitute  the  borohydride  oxidation  mechanism  are  still  a 
matter  of  debate.  Reaction  (R1 )  describes  the  ideal  global  oxidation 
reaction,  but  the  8e~  are  likely  provided  to  the  anode  through  a 
series  of  single  electron  transfers  involving  intermediate  species 
adsorbed  on  the  electrocatalyst  surface  and  in  solution  [11,12],  The 
intermediate  species  provide  opportunities  for  alternative  reaction 
pathways,  some  involving  chemical  steps  which  do  not  yield 
electrons  to  the  anode.  An  analogous  situation  exists  at  the  cathode, 
for  which  Reaction  (R2)  describes  the  ideal  global  H202  reduction 
reaction. 

Despite  a  lack  of  consensus  in  the  literature  as  to  the  detailed 
borohydride  oxidation  mechanism,  a  system  of  three  competing 
reactions  has  been  used  to  describe  the  processes  occurring  at  a 
DBFC  anode.  The  first  reaction  is  (Rl).  The  second  reaction  is 
borohydride  hydrolysis  (Reaction  (R4)),  in  which  the  anode  catalyst 
accelerates  the  chemical  conversion  of  BH4  to  B02 .  Hydrolysis  is 
undesirable  because  it  uselessly  consumes  fuel  and  can  produce 
bubbles  which  inhibit  the  transport  of  BH4  to  the  electrocatalyst 
surface  for  Reaction  (Rl). 

BH4  (aq)  +  2H20— >B02  (aq)  +  4H2  (R4) 

The  third  reaction  is  hydrogen  oxidation  (Reaction  (R5)),  which 
can  consume  H2  produced  by  Reaction  (R4)  and  recover  electrons 
that  would  otherwise  be  lost. 

2  H20  +  2e“  ^  H2  +  20H~  (aq)  fg5  =  -0.828  V  vs.  RHE 

(R5) 

The  net  result  of  Reactions  (Rl),  (R4)  and  (R5)  depends  on  the 
anode  potential  and  species  concentrations  in  solution.  Reaction 
(Rl)  dominates  when  the  BH4  concentration  is  low,  the  hydroxide 
concentration  is  high,  or  the  electrode  potential  is  substantially  less 
negative  than  £°ri  [13],  When  these  conditions  are  not  met, 
adsorbed  hydrogen  can  accumulate  on  the  anode  electrocatalyst, 
shifting  the  net  reaction  in  favor  of  Reaction  (R4)  as  adsorbed 
hydrogen  forms  H2  rather  than  water  [14  .  Reaction  (R5)  contrib¬ 
utes  electrons  to  the  anode  at  potentials  less  negative  than  £°R5  and 
consumes  electrons  at  more  negative  potentials. 

The  undesirable  hydrolysis  reaction  effectively  reduces  the 
number  of  electrons  ne iec  produced  per  BH4  fuel  molecule 
consumed.  Ideally,  the  cell  is  operated  such  that  Reaction  (Rl )  goes 
to  completion  and  ne iec  =  8,  but  many  operating  conditions  lead  to 
intermediates  that  favor  Reaction  (R4)  and  thereby  reduce  ne\ec  and 
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the  power  generated  by  the  DBFC.  The  power  is  reduced  not  only  by 
the  lower  current  but  also  by  the  mixed  potentials  that  arise  due  to 
the  intermediate  species  and  multiple  oxidation  reactions.  Several 
boron-containing  intermediate  species  have  been  proposed.  One 
example  is  BH3OH  ,  whose  oxidation  may  be  rate  limiting  to  Re¬ 
action  (R1 )  at  low  anode  potentials  15]  and  whose  hydrolysis  may 
be  a  key  contributor  to  the  extent  of  Reaction  (R4)  [16], 

The  BH4  oxidation  pathways  are  controlled  in  part  by  the  se¬ 
lection  of  anode  catalyst,  with  Au  and  Pt  having  received  the  most 
attention  [17],  Recent  studies  have  focused  on  alloys  and  nano- 
structured  catalysts  with  improved  activity  and/or  selectivity  for 
BH4  electrochemical  oxidation  [17—19],  The  various  metallic  elec¬ 
trocatalysts  differ  in  their  activity  for  Reactions  (R4)  and  (R5)  with 
Ag  and  Au  being  the  least  active  and  Pt-  and  Pd-based  electro¬ 
catalysts  being  the  most  active  for  both  reactions  1 15,18,20],  The 
relative  rates  for  the  global  Reactions  (Rl)  and  (R4)  with  a  given 
catalyst  depend  on  the  anode  potential,  BH4  concentration  and 
solution  pH  with  increased  pH  improving  selectivity  for  Reaction 
(Rl)  [8],  The  anode  reaction  can  be  approximated  by  the  ideal 
electrochemical  oxidation  (Rl )  for  many  conditions.  For  example, 
Au  anodes,  in  spite  of  their  lower  activity,  are  highly  selective  for 
Reaction  (Rl)  when  the  BH4  concentration  is  low  relative  to  the 
OH~  concentration  [13],  However,  the  same  study  showed  that  BH4 
concentrations  above  20%  of  the  OH~  concentration  lead  to  sig¬ 
nificant  H2  evolution  because  Au  does  not  oxidize  H2  well.  On  more 
active  Pt-  and  Pd-based  electrocatalysts,  the  range  of  operating 
conditions  that  facilitate  the  complete  Reaction  (Rl)  are  limited 
because  these  catalysts  have  higher  activity  than  Au  for  the  hy¬ 
drolysis  reaction.  However,  if  the  cell  is  operated  at  high  current 
densities  where  the  anode  potential  is  less  negative  than  £°R5, 
active  catalysts  such  as  Pt  and  some  of  its  alloys  readily  drive  the 
reverse  of  Reaction  (R5)  [15,21].  Thus,  at  higher  current  densities 
the  full  BH4  oxidation  with  neiec  approaching  eight  can  be  achieved 
[22]  if  the  catalyst  layer  is  designed  to  ensure  good  access  of  OH~  to 
the  catalyst.  High  current  densities  favor  full  oxidation  in  part  by 
maintaining  low  concentrations  of  BH4  near  the  anode  due  to 
transport  limitations  in  the  presence  of  excess  OHA 

The  combination  of  charge-transfer  and  chemical  catalytic  re¬ 
actions  at  the  cathode  can  be  described  by  three  reactions  analo¬ 
gous  to  Reactions  (Rl ),  (R4)  and  (R5)  at  the  anode.  The  first  cathode 
reaction  is  the  complete  reduction  of  H2C>2  in  an  acidic  medium  in 
Reaction  (R2).  The  second  reaction,  (R6),  is  the  heterogeneous 
chemical  decomposition  of  H2O2  into  O2  and  water,  accelerated  by 
the  cathode  electrocatalyst.  Reaction  (R6)  at  the  cathode  is  analo¬ 
gous  to  Reaction  (R4)  at  the  anode. 

2  H202  (aq)^02  +  2H20  (R6) 

Like  Reaction  (R4),  the  H202  decomposition  Reaction  (R6)  re¬ 
duces  the  number  of  electrons  per  reactant  molecule  consumed. 
Similarly,  Reaction  (R6)  can  produce  gaseous  02,  which  can  occlude 
the  cathode  catalyst.  Cathode  design  can  mitigate  this  problem,  and 
at  adequately  low  cathode  potentials,  02  produced  by  Reaction  (R6) 
may  be  reduced  electrochemically  at  the  cathode  by  Reaction  (R7) 
to  recover  the  otherwise  lost  electrons. 

02  +  2H+  +  2e-^2  H202(aq)  ER7  =  0.695  V  vs.  RHE  (R7) 

Experiments  with  DBFC  chemistries  different  from  this  study 
have  been  reported  in  the  literature.  The  fuel  solution  in  a  DBFC 
consists  of  aqueous  BH4  and  OH~  ions  with  an  alkali  metal  cation, 
and  while  Li+  and  I<+  have  been  suggested,  the  most  common  is 
Na+  [9,23],  DBFC  fuel  solutions  include  OH~  as  a  reactant  and  sta¬ 
bilizer  for  BH4,  which  can  be  stored  for  months  in  high  pH  media 
(pH  >13)  [24,25],  A  range  of  oxidizers  has  explored,  including  02  (in 
air  or  pure)  and  H202  [9,10],  When  H202  is  chosen,  it  can  be  stored 


in  acidic  or  alkaline  media.  Acidic  solutions  provide  higher  standard 
half-cell  potential  for  H202  reduction  (1.763  V  vs.  RHE)  than  alka¬ 
line  solutions  (0.870  V  vs.  RHE)  [7],  Alkaline  solutions  offer  ways  to 
mitigate  system-level  problems  such  as  volume  imbalance  due  to 
water  crossover  through  the  membrane  [6],  Among  acidic  media, 
some  anions  impede  Reaction  (R2)  more  than  others,  with  the  best 
activity  occurring  in  the  hydrohalic  acids  HI,  HBr,  and  HC1  [26,27], 
However,  a  common  choice  among  experiments  reported  in  the 
literature  is  H2SO4. 

Alternative  cell  geometries  and  topologies  have  also  been 
explored.  Experimental  results  have  been  published  for  cells  with  a 
catalyst  layer  on  the  channel  wall  [1,6],  catalyst  coated  on  the 
membrane  [23,28,29  and  catalyst  distributed  throughout  the 
channel  on  a  porous  medium  [4,30],  The  relative  strengths  and 
weaknesses  of  the  different  electrocatalyst  locations  for  both  the 
anode  and  the  cathode  are  not  clear  because  of  the  lack  of  good 
models  or  direct  experimental  comparison. 

While  there  has  been  significant  work  on  the  electrochemistry 
of  DBFCs,  practical  DBFC  cell  design  and  operating  spaces  have  not 
been  explored  as  thoroughly.  Experiments  on  fundamental  kinetics 
as  well  as  electrode  geometry  have  provided  some  insight  into 
what  is  feasible  for  DBFCs,  but  relationships  between  cell  designs 
and  performance  are  not  understood.  The  complex  coupling  of  the 
reaction  pathways  with  transport  through  reactant  concentrations, 
particularly  for  the  anode,  raise  many  questions  as  to  how  to 
translate  the  numerous  electrocatalyst  kinetic  studies  into  pre¬ 
dicting  cell  performance  for  practical  design.  The  range  of  reactant 
and  product  concentrations  along  the  flow  path  of  a  practical  cell 
with  substantial  reactant  utilization  can  have  a  significant  impact 
on  reaction  overpotentials  and  transport  limitations,  and  thus 
overall  cell  efficiencies  and  power  densities.  Understanding  such 
relationships  is  critical  for  practical  DBFCs  to  achieve  their  prom¬ 
ising  potential.  In  this  regard,  carefully  constructed  numerical 
models  of  DBFCs  can  explore  cell  configurations  and  designs  that 
provide  preferred  conditions  for  coupling  complex  transport  with 
the  kinetics  to  favor  the  highest  power  densities  and  the  most 
effective  conversion  of  fuel  and  oxidizer  to  useful  energy  (i.e.,  ef¬ 
ficiency).  A  model  can  also  provide  state  information  that  is  useful 
for  understanding  how  the  cell  operates,  but  is  otherwise  difficult 
to  measure.  This  study  presents  a  2-D  steady  state  finite  volume 
model  to  explore  how  cell  geometry  and  operating  conditions 
impact  DBFC  performance  under  conditions  which  favor  Reactions 
(Rl)  and  (R2).  The  specific  goals  of  this  work  were: 

1.  Explore  the  performance  that  may  become  available  if  efforts  to 

improve  DBFC  electrocatalysts  and  membranes  are  successful. 

2.  Understand  trends  relating  DBFC  cell  design  parameters  to 

performance,  so  that  cells  can  be  designed  to  favor  conditions  in 

which  Reactions  (Rl)  and  (R2)  predominate. 

The  modeled  system  is  ideal  in  that  only  Reactions  (Rl )  and  (R2) 
take  place,  and  the  membrane  is  permeable  only  to  Na+  and  water. 
Realistically,  improved  electrocatalysts  and  membranes  are  un¬ 
likely  to  eliminate  Reactions  (R4)— (R7)  and  the  transport  of  species 
other  than  Na+  through  the  membrane.  The  results  presented  here 
represent  the  performance  of  an  ideal  DBFC,  which  realistic  DBFC 
performance  may  approach  as  the  technology  improves. 

DBFC  models  have  been  published  prior  to  this  work.  Verma  and 
Basu  [31  ]  modeled  an  02  DBFC  consisting  of  a  Nafion  membrane  in 
contact  with  porous  catalyst  layers  and  a  well-mixed  volume  of 
fuel.  Mass  transport  to  the  cathode  was  governed  by  Fields  law  (02 
diffusion  in  air)  and  transport  to  the  anode  was  neglected.  All 
charge  transport  processes  in  the  cell  were  described  by  a  fit  to 
measurements  of  ohmic  resistance.  Sanli  et  al.  [32]  published  a 
similar  model,  but  with  H202  as  the  oxidizer  and  mass  transport 
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neglected  entirely.  Both  [31]  and  [32]  ignore  down-the-channel 
effects.  Shah  et  al.  [33]  published  the  most  detailed  DBFC  model 
to  date,  for  a  cell  topology  similar  to  that  of  PEMFCs.  The  cell  con¬ 
sisted  of  a  Nafion  membrane  with  a  porous  catalyst  layer  on  each 
side;  the  catalyst  layers  were  separated  from  reactant  flow  channels 
by  porous  diffusion  layers.  The  model  predicts  cell  voltage  at  a 
given  current  density  by  calculating  activation  and  ohmic  over¬ 
potentials  explicitly  as  functions  of  current  density,  and  then  sub¬ 
tracting  the  overpotentials  from  the  open  circuit  voltage. 
Concentration  overpotentials  were  included  by  estimating  con¬ 
centrations  near  the  electrodes,  and  using  those  concentrations  in 
the  calculation  of  activation  overpotentials.  Species  concentrations 
were  assumed  to  be  uniform  in  the  flow  channels,  but  transport 
from  the  flow  channels  to  the  catalyst  layer  was  approximated  by  a 
Nernst  diffusion  layer  thickness  correlation.  Down-the-channel 
effects  were  not  addressed.  Of  prior  DBFC  models,  only  one  [34] 
included  transport  in  a  meaningful  way  and  none  addressed 
down-the-channel  effects  on  cell  performance.  No  prior  DBFC 
model  has  studied  the  geometry  depicted  in  Fig.  1. 

Other  fuel  cell  topologies  and  chemistries  have  been  studied 
using  2-D  finite  volume  models.  A  relevant  example  is  that  of 
Sprague  and  Dutta  [35,36],  who  modeled  a  fuel  cell  similar  to  Fig.  1 
with  the  exception  that  it  had  no  membrane.  The  cell  had  a  single 
liquid  electrolyte  and  electrically  neutral  fuel  and  oxidizer  species 
which  reacted  in  a  generalized  way  governed  by  the  Frumkin— 
Butler— Volmer  equation.  The  Poisson— Nernst— Plank  equations 
were  solved  to  resolve  electrochemical  double  layers  near  the 
electrodes  and  the  distributions  of  species  and  charge  in  the 
channel.  The  model  presented  here  is  solved  in  a  similar  way,  but  it 
differs  from  that  of  Sprague  and  Dutta  in  several  key  aspects.  The 
model  in  this  study  includes  a  membrane,  has  different  electrolytes 
and  chemistries  (specific  to  a  DBFC)  at  the  anode  and  cathode,  and 
solves  for  the  electric  field  by  explicitly  enforcing  electroneutrality 
rather  than  by  solving  Poisson’s  electrostatic  equation. 

This  work  goes  beyond  previous  DBFC  modeling  studies  by 
examining  a  new  cell  topology,  including  down-the-channel  effects 
with  a  2-D  model,  and  exploring  the  impacts  of  cell  geometry  and 
operating  parameters  on  performance.  The  model  here  provides  a 
framework  in  which  to  study  cell  performance  under  conditions 
where  complete  borohydride  oxidation  and  complete  hydrogen 
peroxide  reduction  are  good  approximations  for  the  electrode  re¬ 
actions.  The  model  can  provide  DBFC  developers  with  guidance  on 
how  the  complex  coupling  of  transport  and  chemistry  impacts 
practical  cell  design.  This  study  establishes  a  basis  for  further 
studies  incorporating  more  realistic  chemistry  and  alternative  cell 
configurations  that  can  accelerate  practical  DBFC  development. 

2.  Mathematical  model 

The  steady  state  performance  of  a  DBFC  is  modeled  with  a  2-D 
finite  volume  approach.  The  model  domain  includes  two  rectan¬ 
gular  channels  separated  by  a  membrane  and  channel  walls  which 
function  as  electrodes  (see  Fig.  1 ).  An  aqueous  fuel  solution  containing 
Na+,  OH~,  BH4 ,  and  BO2  enters  the  anode  channel  and  an  aqueous 
oxidizer  solution  containing  H2O2,  H+,  and  SC>4~  enters  the  cathode 
channel.  Both  solutions  flow  in  the  same  direction.  Charge  transfer 
reactions  take  place  at  the  channel  walls  while  Na+  passes  through  a 
Nafion  membrane  in  the  Na+  form  from  fuel  channel  to  oxidizer 
channel. 

The  operating  conditions  and  cell  geometry  in  this  study  enable 
several  simplifications.  The  fuel  and  oxidizer  solutions  are  liquids 
and  therefore  effectively  incompressible.  The  flows  in  the  reactant 
channels  are  laminar  because  the  Reynolds  numbers  are  less  than 
150.  The  fuel  and  oxidizer  solutions  are  electrically  neutral  because 
the  electrochemical  double  layers  are  thin  and  treated  as  part  of  the 


electrode  interfaces.  The  electrodes  are  good  electrical  conductors 
and  therefore  all  points  along  the  channel  for  each  electrode  share 
the  same  electric  potential.  Homogeneous  reactions  may  take  place 
in  the  fuel  and  oxidizer  solutions,  but  at  high  pH  (pH  >13)  in  the  fuel 
and  low  pH  (pH  <1)  in  the  oxidizer  they  occur  so  slowly  in  com¬ 
parison  to  the  electrode  reactions  that  they  can  be  neglected  without 
significantly  impacting  the  results  [24,25],  The  membrane  is  fully 
hydrated  Nafion  115  (145  pm  thick)  in  the  Na+  form,  which  is  thin 
enough  to  make  x-direction  fluxes  in  the  membrane  insignificant. 

Four  simplifying  assumptions  accommodate  unknown  trans¬ 
port  parameters: 

•  All  activity  coefficients  are  1,  i.e.  the  fuel  and  oxidizer  solutions 
are  ideal  electrolytes. 

•  Only  Na+  and  H20  pass  through  the  membrane  (negligible  BH4 
crossover). 

•  Solution  viscosity  is  the  same  as  pure  water. 

•  Heat  is  removed  quickly  enough  to  maintain  isothermal 
conditions. 


2.1.  Model  geometry,  state  variables  and  solution  approach 

The  model  domain  is  subdivided  into  channels  and  interfaces. 
Channels  are  discretized  into  numerical  cells  with  state  variables 
for  pressure,  x-  and  y-direction  mass-averaged  velocities,  electric 
potential,  and  species  mass  fraction.  Nomenclature  for  the  state 
variables  and  other  model  parameters  is  provided  in  Appendix  A. 
Fig.  1  shows  orientations  of  the  x-  and  y-coordinates.  Velocity 
storage  locations  are  staggered  from  scalar  variables  by  one-half 
discretization  to  avoid  odd— even  decoupling.  All  interfaces  are 
broken  into  discrete  line  segments  (numerical  cells)  having  state 
variables  for  electric  potential  and  species  mass  fractions.  The 
membrane  interface  also  stores  pressure. 

The  numerical  grid  density  varies  linearly  in  the  x-  and  y-di- 
rections  with  the  highest  density  near  the  inlet,  electrodes  and 
membranes  to  capture  the  development  of  entrance  effects  and 
boundary  layers.  Trial  and  error  showed  that  consistent  results  are 
achieved  when  the  smallest  x-discretization  is  0.25  times  the  mean 
x-discretization  length  and  the  smallest  y-discretization  is  0.10  times 
the  mean  y-discretization  length.  The  grid  density  was  increased 
until  the  solution  achieved  grid  independence.  Results  presented 
here  require  49,800  state  variables  to  simulate  a  200  mm  long 
channel  broken  into  150  x-discretizations  and  17  y-discretizations. 

A  solution  is  the  set  of  state  variable  values  which  conserve 
mass,  momentum,  species  and  charge  and  are  consistent  with  the 
boundary  conditions.  From  a  top-level  point  of  view,  the  total  cell 
current  is  calculated  for  specified  cell  voltage,  geometry  and  oper¬ 
ating  conditions.  The  model  is  set  up  in  MATLAB  and  solved  using  a 
Newton  solver  in  K1NSOL.  KINSOL  is  a  non-linear  equation  solver 
available  from  the  Lawrence  Livermore  National  Laboratory  as  part 
of  a  collection  of  equation  solving  tools  called  SUNDIALS  [37],  The 
MATLAB  code  sets  up  the  model  geometry,  discretizes  the  domain, 
applies  boundary  conditions,  and  provides  a  function  KINSOL  calls 
to  evaluate  the  residuals  (errors)  associated  with  each  prospective 
solution.  KINSOL  runs  until  the  residual  2-norm  is  driven  below  a 
specified  threshold  (typically  1CT3),  at  which  point  the  solution  has 
been  found.  Solution  time  using  one  3.3  GHz  Xeon  processor  core  is 
typically  ~80  min,  but  it  is  reduced  to  ~40  min  when  the  initial 
guess  is  from  a  solution  for  similar  conditions. 

2.2.  Governing  equations  and  boundary  conditions 

The  residual  for  each  state  variable  is  the  error  in  the  associated 
governing  equation.  The  state  variables  require  governing 
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equations:  conservation  of  mass  (associated  with  pressure),  con¬ 
servation  of  momentum  (associated  with  velocity),  conservation  of 
species  (associated  with  mass  fraction)  and  electroneutrality 
(associated  with  electric  potential).  Boundary  conditions  are 
applied  at  the  inlets,  outlets,  electrodes,  electrode  interfaces  and 
membrane  interfaces.  Some  variables  have  one  simple  boundary 
condition  and  others  have  multiple  mixed  boundary  conditions 
which  are  linked  across  interfaces  by  the  related  fluxes. 

2.2. 2.  Pressure 

Pressure  is  associated  with  mass  conservation  because  mass 
fluxes  are  predominantly  driven  by  pressure  gradients.  At  steady 
state,  mass  conservation  in  the  channels  is  the  sum  of  mass  fluxes 
across  the  boundaries  of  a  numerical  cell,  therefore  the  residuals 
associated  with  mass  conservation  are: 

Kh  =  -v-T  (i) 

Mass  conservation  is  also  used  to  solve  for  pressure  at  the 
membrane  interfaces,  where  two  fluxes  in  the  y-direction  are 
involved  (one  mass  flux  from  the  channel  and  one  from  the 
membrane): 


Table  1 

Velocity  boundary  conditions.  The  x-coordinate  runs  from  inlet  to  outlet  and  the  y- 
coordinate  runs  from  electrode  to  membrane,  as  in  Fig.  1. 


Boundary 

Description 

Boundary  condition(s) 

Channel 

Fully  developed 

vx  =  1.5  V(1  -  (2y//i  -  l)2) 

inlets 

with  specified 
volumetric  flow  rate 

Vy  =  0 

Channel 

outlets 

Fully  developed 

dvjdx  =  0,  dVy/dX  =  0, 

Electrode 

No  slip  condition 

vx  =  0 

interface 

Impermeable  boundary 

Vy  =  0 

Membrane 

No  slip  condition 

vx  =  0 

interface 

No  gradient  in  membrane 

dVy/dy  =  0 

dictated  by  the  respective  flux  matches.  Ultimately,  the  mass  frac¬ 
tions  on  each  side  of  the  membrane  are  related  by  the  membrane 
species  fluxes,  and  mass  fractions  at  the  electrodes  are  related  to 
the  electrode  potential  boundary  conditions  by  the  reaction  rates. 

2.2.3.  Velocity 

Momentum  conservation  in  the  channels  can  be  written  as  the 
sum  of  advection,  pressure,  diffusion  and  electric  body  force  terms: 


D  P 

iVm_int 


Jch-  J 


(2) 


The  x-direction  pressure  boundary  condition  is  constant  pres¬ 
sure  at  the  outlets  (P  =  100  kPa).  The  y-direction  pressure  boundary 
conditions  are  zero  pressure  gradient  at  the  electrode  (0P/3y  =  0)  in 
recognition  of  the  zero  net  mass  flux  there. 


2.2.2.  Mass  fractions 

Species  concentrations  are  derived  from  species  conservation, 
which  is  written  in  terms  of  mass  fraction  to  simplify  the  mass  and 
momentum  conservation  equations.  The  concentration  of  species  k 
is  calculated  from  mass  fraction  by  Q  =  pYjJW^.  The  species  con¬ 
servation  equation  is  the  sum  of  species  flux  and  storage  terms: 
0  =  V  •  J  k  +  3 (pYk) /dt.  The  storage  term  can  be  rearranged  to 
isolate  dYiJdt  on  the  left  hand  side,  which  should  be  zero  at  steady 
state  and  becomes  the  mass  fraction  residual.  The  mass  storage 
term  Pch  =  3p/3t  is  zero  at  steady  state,  but  it  is  retained  to  improve 
solver  convergence  by  maintaining  the  link  between  species  and 
mass  conservation.  The  mass  fraction  residual  for  species  k  in  the 
channel  is  therefore: 

Rl?ch=-p(-V-7k-YkRPch)  (3) 

Species  conservation  is  also  used  to  find  mass  fraction  residuals 
at  the  electrode  and  membrane  interfaces,  but  it  becomes  a  flux 
match  because  the  interfaces  have  only  two  species  fluxes  and  no 
mass  storage  term.  At  the  electrode  interfaces  there  is  one  flux  from 
the  channel  and  one  flux  from  the  electrode  (due  to  reactions  there) 
for  each  species  k: 

<L  =  Ik.ch  -  lk,m  (4) 

Mass  fraction  residuals  at  the  membrane  interfaces  involve  one 
flux  from  the  channel  and  one  through  the  membrane  for  each 
species  k. 


3(p  v  )  =  _ Vplf .  q?  _  vP  +  /xV2T?  -  pcV0  (6) 

The  electric  body  force  term  pcV0  accounts  for  the  influence  of 
the  electric  potential  gradient  on  ions  in  solution.  Taking  only 
components  involved  in  the  x-direction  momentum  balance  and 
rearranging  to  obtain  3vx/3t  (which  should  be  zero  at  steady  state) 
yields  the  x-direction  velocity  residual  in  the  channels: 


pxM 

Kch 


1 

P 


3 

3x 


(PFx|i'x|) 


3 y 


(pVXVy) 


3  P 
3x 


+  P 


(7) 


The  mass  conservation  storage  term  Rch  is  retained  in  equation 
(7)  for  the  same  reasons  as  in  equation  (3).  The  y-direction  velocity 
residuals  R have  the  same  form  as  equation  (7): 


oyM  _  1 
ch  -  7 
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(8) 

Velocity  and  pressure  gradients  in  the  x-direction  are  estimated 
by  upwind  differencing  because  transport  in  the  x-direction  is 
dominated  by  advection,  and  upstream  state  variables  have  more 
influence  over  advection  than  downstream  state  variables  [38], 
Velocity  and  pressure  gradients  in  the  y-direction  are  estimated  by 
center  differencing  because  they  are  small  and  transport  in  the  y- 
direction  is  dominated  by  diffusion  and  migration.  Velocity 
boundary  conditions  were  assigned  at  the  channel  inlets,  outlets 
and  interfaces  as  shown  in  Table  1. 


R 


MF 
k.m  i 


J  k.ch  J  i 


(5) 


The  x-direction  mass  fraction  boundary  conditions  are  specified 
at  the  inlet  for  each  species  k  based  on  the  fuel  and  oxidizer  con¬ 
centrations.  The  y-direction  mass  fraction  boundary  conditions  are 
the  mass  fractions  at  the  electrode  and  membrane  interfaces 


2.2.4.  Electric  potential 

The  electric  potential  in  an  electrolyte  solution  is  related  to  the 
net  charge  density  by  Poisson’s  electrostatic  equation. 

V2d>  =  — ^  (9) 
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The  small  permittivity  of  free  space  (eo  =  8.85419  x  1CT12  F  rrT1 ) 
makes  the  electric  field  gradient  sensitive  to  any  deviation  from 
electroneutrality,  and  in  real  systems  such  deviations  only  occur  in 
electrochemical  double  layers.  Elsewhere  the  net  charge  density 
can  be  assumed  to  be  zero  [39],  Some  models  [35]  have  solved 
equation  (9)  in  order  to  resolve  double  layers,  but  in  this  work  the 
double  layers  are  considered  part  of  the  interfaces,  so  the  electric 
potential  is  found  by  enforcing  electroneutrality.  The  electric  po¬ 
tential  residual  is  equal  to  the  net  charge  density  (i.e.  deviation 
from  electroneutrality)  and  is  given  by  equation  (10),  which  is 
solved  to  find  the  electric  potential  in  the  channels  and  at  the 
interfaces. 


Table  3 

Thermophysical  and  transport  properties  in  fully  hydrated  Nafion  115  in  the  Na+ 
form  at  298  K. 


Property 

Value 

Na+  mobility,  uNa+  mpm 

2.7  X  to-8  m2  V-1  s-1 

[52] 

H2O  permeability,  pHz 0 

1.7  x  1CT14  m  Pa-1  s“ 

‘[53] 

Na+  concentration, 

^Na+  =  ^SOj 

1.13  kmol  itt3  [54] 

Na+— H20  electro-osmotic 
drag  coefficient  in  fully 
hydrated  Nafion  in  the 

Na+  form,  nd 

9.2  kmol  H20/kmol  Na 

+  [52] 

&  =  Fj2zkCk  (10) 

k 

Charge  conservation  implies  that  the  total  anode  current  must 
equal  the  total  cathode  current  at  steady  state.  The  electro¬ 
neutrality  condition  ensures  the  total  electrode  currents  are  equal 
by  conserving  charge  throughout  the  model  domain. 

The  x-direction  boundary  conditions  on  electric  potential  are 
zero  electric  field  at  the  inlets  and  outlets  (di f>/3x  =  0)  and  the  y- 
direction  boundary  conditions  are  the  electrode  potentials.  The 
anode  electric  potential  is  zero  to  establish  a  reference  for  the 
system,  and  the  cathode  electric  potential  is  the  specified  cell 
voltage. 

2.3.  Species,  mass  and  charge  fluxes  in  the  channels 

The  species,  mass  and  charge  fluxes  crossing  scalar  numerical 
cell  boundaries  must  be  calculated  in  order  to  evaluate  the  related 
conservation  equation  residuals.  Solute  species  mole  fluxes  in  the 
channel  are  described  by  the  Nernst— Plank  equation  (equation 
(11))  which  includes  contributions  from  diffusion,  migration  and 
advection. 

Jk  =  -DkVCk  -  zkukFCkV(j)  +  pYk~v /Wk  (11) 


The  total  charge  flux  is  a  charge-weighted  sum  of  the  mole 
fluxes: 

Tc  =  fJ2zi<7i<  (14) 

k 

Equations  (9)— (14)  require  concentration  and  electric  potential 
gradients  at  the  discretization  boundaries.  Gradients  are  estimated 
by  center  differencing  because  it  tends  to  be  more  numerically 
stable  for  diffusion  and  migration  [38],  Table  2  lists  the  transport 
and  thermophysical  properties  appearing  in  equations  (11)— (14). 

2.4.  Species  and  charge  fluxes  through  the  membrane 

A  1-D  sub-model  was  used  to  calculate  the  Na+  and  H2O  mole 
fluxes  across  a  fully  hydrated  Nafion  115  membrane  at  each  point 
down  the  channel.  BH2  crossover  through  Nafion  membranes  has 
been  observed  in  experiments  [5,40],  but  was  omitted  because  BH4 
migration  fluxes  are  oriented  from  cathode  to  anode,  so  the  rate  of 
crossover  should  be  small  at  operating  conditions  far  from  OCV  [40] 
where  the  electric  potential  gradient  in  the  membrane  is  large. 
Table  3  lists  transport  parameters  in  the  membrane  sub-model. 

Diffusion  and  migration  drive  Na+  through  the  membrane,  so 
the  net  Na+  mole  flux  from  fuel  channel  to  oxidizer  channel  is  given 
by  equation  (15): 


Since  diffusion  cannot  produce  a  net  mass  flux,  the  water 
diffusion  mass  flux  is  made  equal  to  the  negative  sum  of  solute 
diffusion  mass  fluxes.  Water  has  no  net  charge,  so  the  migration 
term  is  dropped  from  the  total  water  mole  flux  in  equation  (12). 

(12) 

The  total  mass  flux  is  the  molecular  mass  weighted  sum  of  the 
mole  fluxes: 

7  =  J2w'<7k  (i3) 
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J  H20  — 


W, 


h2o 


E  (Wk  +  pVW 

fc^H20 


Jl 


=  -D» 


dCNa  d(f> 

a~dy  zNaUNaK-Na-j^ 


(15) 


The  concentration  and  electric  potential  gradients  are  found  by 
assuming  linear  profiles  in  the  y-direction  between  the  membrane 
interfaces.  To  maintain  electroneutrality,  the  concentration  of  Na+ 
must  be  equal  to  the  concentration  of  SO3  groups:  CNa+  =  Csor- 
Water  can  cross  the  membrane  due  to  diffusion,  permeation  and 
electro-osmotic  drag  induced  by  the  Na+  flux.  Water  diffusion  is 
neglected  because  the  concentrations  on  each  side  of  the  mem¬ 
brane  are  nearly  the  same.  Thus,  the  net  water  flux  is  determined  by 
the  permeation  and  electro-osmotic  drag  terms. 


J  h2o.i 


Ph2 ii7  +  nd  J  Na.mem 

wH20 


(16) 


Table  2 

Diffusivities  and  apparent  molar  volumes  of  fuel  and  oxidizer  solutes. 


Species 

Diffusivity  in  H20  at 
infinite  dilution  and 

298  I<  [m2  s-’  x  109] 

Apparent  volume  in 

H20  at  infinite  dilution 

and  298  K 

[m3  kmol-1  x  103] 

BH4 

2.42  [46] 

— 5.83a 

bo2 

0.814  [47] 

-14.5  [48] 

H+ 

9.312  [39] 

0  [48] 

H202 

1.19  [49] 

22.17  [49] 

Na+ 

1.334  [39] 

-1.11  [48] 

OH" 

5.260  [39] 

-4.18  [48] 

sol- 

0.625  [50] 

24.8  [51  ] 

1  Value  unavailable  but  estimated  to  be  smaller  than  B02  by  the  diffusivity  ratio. 


2.5.  Species  and  charge  fluxes  at  the  electrodes 

Multiple  charge  transfer  reactions  can  take  place  at  an  electrode 
simultaneously;  the  rate  of  each  reaction  can  be  modeled  as  the 
sum  of  an  anodic  reaction  rate  supplying  electrons  to  the  electrode 
and  a  cathodic  rate  removing  electrons.  The  rate  r  of  reaction  q  is 
given  by  equation  (17),  in  which  subscripts  a  and  c  refer  to  anodic 
and  cathodic  rate  parameters,  respectively. 

rq  =  (17) 

k  k 
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Both  terms  in  the  rate  expression  depend  on  a  rate  constant  k, 
local  species  concentrations  C/(,  reaction  rate  orders  with  respect  to 
each  species  g/(,  number  of  electrons  transferred  in  the  rate  limiting 
step  ne,  charge  transfer  symmetry  factor  /?  and  the  electric  potential 
difference  between  the  electrode  and  electrolyte  A <p. 

The  net  flux  of  species  k  at  the  electrode  is  needed  to  solve  the 
species  conservation  equation  at  the  electrode  interface  (equation 
(4)).  The  net  species  flux  is  the  sum  of  contributions  from  each 
reaction  occurring  at  the  electrode  as  shown  in  equation  (18);  each 
contribution  depends  on  a  reaction  rate  rq,  stoichiometry  v^q  and 
roughness  factor  5.  The  roughness  factor  accounts  for  the  difference 
between  electrochemical  and  geometric  electrode  surface  area. 

J  fc,rxn  =  ^  y  ]  rQvk.a  (IS) 

4 

The  net  current  density  at  the  electrode  is  found  by  summing 
the  current  density  contributions  from  each  species: 

!net  —  F  y  ]  J  kjxnZk  (19) 

k 

An  advantage  of  this  approach  is  that  it  accommodates  mixed 
potentials,  which  is  not  possible  with  the  traditional  Butler— 
Volmer  approach  because  the  equilibrium  potential  cannot  be 
predicted  by  thermodynamics.  The  present  work  assumes  only  one 
reaction  takes  place  at  each  electrode,  but  the  model  could  be 
expanded  to  describe  the  influence  of  competing  reactions  and 
intermediates  if  the  rate  parameters  were  known. 

A  disadvantage  of  this  approach  is  that  it  requires  rate  param¬ 
eters  (k  and  (1),  which  are  rarely  known.  Since  the  goal  of  this  work 
is  to  identify  performance  trends,  the  reaction  rate  values  are  less 
important  than  capturing  the  correct  functional  dependencies.  The 
following  educated  guesses  were  made  for  reaction  rate 
parameters: 

•  neq  is  set  to  1  because  it  is  assumed  the  overall  rate  is  deter¬ 
mined  by  an  intermediate  step  which  transfers  one  electron. 

•  1 8q<a  and  (lqx  are  set  to  0.5  —  a  common  assumption  in  the 
absence  of  known  values.  [41  ]. 

•  For  the  borohydride  oxidation  reaction  (Reaction  (Rl)),  the 
anodic  rate  is  assumed  to  be  first  order  with  respect  to  Cbh:  and 
Com  and  the  cathodic  rate  is  assumed  to  be  first  order  with 
respect  to  Cbo2 .  There  is  experimental  evidence  for  the  first  or¬ 
der  dependency  on  Cbht  [42], 

•  For  the  hydrogen  peroxide  reduction  reaction  (Reaction  (R2)), 
the  anodic  rate  is  assumed  to  be  first  order  with  respect  to  CH+ 
and  CH2o2  and  the  cathodic  rate  has  no  concentration 
dependence. 


2.6.  Composition  equation  of  state 

The  fuel  and  oxidizer  solutions  are  incompressible,  but  an 
equation  of  state  is  still  needed  to  relate  solution  mass  density  to 
composition.  A  state  equation  (equation  (20))  relates  solution  mass 


Table  4 

Baseline  case  and  variations  from  the  baseline. 


Parameter 

Baseline 

Variations  from  baseline 

Anode  inlet  [NaBH4] 

0.3  M 

0.1,  0.2,  0.4,  0.5  M 

Channel  depth 

1  mm 

0.50,  0.75  mm 

Membrane  thickness 

145  pm 

1.45,  72.50  pm 

Inlet  fuel  flow  rate 

60  mL  min-1 

15,  30  mL  min-1 

Forward  reaction 

rate  constants 

106  m4  mol-1  s_1 

10°,  103  m4  mol-1  s_1 

density  to  mass  fraction  by  accounting  for  the  apparent  molar 
volume  of  each  solute. 

_  _ Ph2o _  f2ni 

P  ~  YH2o  +  pH20j:k(5kYk/Wk)  1  J 

Pn2o  varies  with  temperature  via  a  polynomial  fit  to  empirical 
data  at  100  kPa. 

3.  Results  and  discussion 

This  study  consists  of  a  baseline  case  and  12  alternative  cases  in 
which  one  parameter  is  altered  (see  Table  4)  to  identify  the  effects  it 
has  on  DBFC  performance.  Many  parameters  remained  the  same  in 
all  cases  (see  Table  5).  The  model  was  run  repeatedly  for  each  case, 
first  with  a  cell  voltage  of  3.4  V  and  then  decreasing  the  cell  voltage 
in  0.1  V  steps  until  it  reached  the  transport  limited  current  density. 
Each  run  yielded  values  for  all  state  variables  in  the  model  domain 
and  post  processing  provided  the  mass,  momentum  and  charge 
fluxes  (i.e.  current  density).  The  state  variables  and  fluxes  were 
used  to  estimate  performance  metrics  such  as  voltage  efficiency, 
power  density  and  fuel  utilization. 

Baseline  values  for  the  channel  dimensions,  membrane  thick¬ 
ness  and  inlet  flow  rates  were  based  on  typical  values  appearing  in 
published  experimental  studies.  The  inlet  BH4  concentrations  were 
chosen  to  operate  in  a  regime  where  complete  borohydride 
oxidation  (Reaction  (Rl))  is  likely  to  dominate  the  anode  reactions 
(at  cell  voltages  other  than  open  circuit).  Fast  forward  reaction  rate 
constants  (klA  at  the  anode  and  k i,c  at  the  cathode)  were  chosen  to 
model  the  performance  of  a  cell  with  an  advanced  electrocatalyst 
and  emphasize  the  influence  of  transport. 

3 A.  Baseline  case 

Voltage  efficiency  is  the  ratio  of  actual  cell  voltage  to  the  ideal 
open  circuit  voltage  (OCV)  predicted  by  thermodynamics.  The 
voltage  losses  in  a  DBFC  stem  primarily  from  the  following:  1) 
competing  electrochemical  reactions  that  depress  OCV  by  shifting 
the  equilibrium  potentials  of  the  anode  and  cathode,  2)  activation 
overpotentials  to  drive  electrons  across  electrode-solution  in¬ 
terfaces,  3)  changes  in  concentration  near  the  electrodes  as  con¬ 
centration  gradients  develop  to  drive  fluxes  to  or  from  the 
electrodes,  and  4)  ohmic  overpotentials  to  drive  charge  from  anode 
to  cathode. 

The  voltage  losses  (as  in  all  fuel  cells)  increase  with  cell  current 
density,  which  can  be  visualized  by  a  polarization  curve;  the 
calculated  baseline  polarization  curve  is  shown  in  Fig.  2(a).  The  cell 
voltages  in  Fig.  2(a)  are  useful  in  a  relative  sense  because  the  model 
neglects  fuel  crossover  and  competing  electrochemical  reactions, 
but  the  slopes  of  the  polarization  curve  in  the  linear  region  at  in¬ 
termediate  current  densities  and  in  the  transport-limited  region  at 
the  highest  current  densities  are  credible  because  they  are  largely 
dictated  by  known  transport  parameters.  The  product  of  cell 
voltage  and  average  current  density  gives  an  average  cell  area  po¬ 
wer  density  (W  itT2  of  electrolyte  membrane),  which  provides  a 


Table  5 

Parameters  common  to  all  cases. 


Parameter 

Value 

Anode  inlet  [NaOH] 

4.0  M 

Anode  inlet  [NaB02] 

10-6  M 

Channel  width 

10  mm 

Temperature 

298  K 

Cathode  inlet  [H2SO4] 

4.0  M 

Cathode  inlet  [H202] 

4.0  M 

Oxidizer  inlet  flow  rate 

60  mL  min 
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Fig.  2.  (a)  Baseline  polarization  curve  showing  activation,  ohmic  and  concentration  overpotentials,  (b)  Baseline  power  density  vs.  cell  voltage.  Points  to  the  right  of  2.2  V  offer  the 
best  combination  of  power  density  and  voltage  efficiency,  (c)  Area  specific  resistance  (ASR)  associated  with  each  overpotential  at  points  down  the  channel,  in  the  baseline  case  at 
2.2  V.  (d)  Baseline  current  density  at  2.2  V  at  points  down  the  channel,  which  correlates  with  the  rate  of  reactant  consumption.  Anode  and  cathode  current  densities  at  each  point 
are  the  same  to  within  0.01%. 


link  to  stack  (volume)  power  density  (W  uT3).  Plotting  cell  area 
power  density  vs.  cell  voltage  as  in  Fig.  2(b)  illustrates  the  tradeoffs 
between  power  density  and  voltage  efficiency.  Fig.  2(b)  shows 
baseline  peak  power  density  occurs  at  2.2  V.  Operating  at  cell 
voltages  higher  than  2.2  V  offers  the  best  combination  of  power 
density  and  efficiency. 

The  voltage  losses  due  to  activation,  concentration  and  ohmic 
resistance  vary  along  the  channel  with  local  concentrations  and 
current  density.  The  magnitude  of  each  loss  can  be  characterized  by 
the  associated  area  specific  resistance  (ASR).  Fig.  2(c)  plots  the 
variation  of  the  ASRs  for  each  electrode  and  the  membrane  along 
the  channel  at  a  cell  voltage  of  2.2  V.  For  the  baseline  reaction  rates, 
the  ASRs  for  concentration  in  the  fuel  channel  and  cathode  acti¬ 
vation  are  largest,  indicating  they  contribute  most  to  losses  in  the 
fuel  cell  at  this  operating  point.  Both  activation  and  concentration 
ASRs  increase  down  the  channel  as  concentration  boundary  layers 
lower  the  availability  of  reactants  at  the  electrodes  and  depress  the 
actual  OCV  further  from  the  ideal  OCV.  Fig.  2(c)  shows  the  impor¬ 
tance  of  variation  down  the  channel  in  a  DBFC,  which  can  only  be 
captured  with  a  multi-dimensional  model. 

The  variation  down  the  channel  in  the  anode  concentration  ASR 
is  explained  in  Fig.  2(d),  which  shows  current  density  as  a  function 


of  distance  along  the  channel.  BH4  is  consumed  rapidly  near  the 
channel  inlet  but  the  consumption  rate  falls  at  points  approaching 
the  outlet.  Thus,  while  channel  fuel  utilization  (fuel  consumed  over 
fuel  in)  increases  with  lengthening  the  channels,  the  increased 
length  results  in  lower  power  densities.  The  channel  fuel  utilization 
in  the  case  shown  in  Fig.  2(d)  is  4.4%.  As  seen  in  Fig.  2(c),  the  fuel 
channel  concentration  ASR  becomes  higher  than  other  resistances 
with  distance  along  the  channel  due  to  the  growing  thickness  in  the 
fuel  concentration  boundary  layer.  Transport  through  this  bound¬ 
ary  layer  in  the  y-direction  is  complicated  because  advection, 
diffusion  and  migration  all  generate  BH4  fluxes  in  the  y-direction. 
Fig.  3  shows  anode  channel  BH4  fluxes  in  the  y-direction,  mid-way 
between  the  inlet  and  outlet,  at  a  cell  voltage  of  2.2  V.  Transport  at 
the  anode  surface  is  dominated  by  diffusion  because  the  concen¬ 
tration  of  BHT  there  is  nearly  zero,  but  migration  plays  a  significant 
role  in  the  center  of  the  channel  where  the  concentration  is  higher 
and  the  electric  field  drives  BH4  toward  the  anode. 

Transport  of  other  species  also  impacts  the  DBFC  performance. 
In  particular,  transport  of  Na+  and  H2O  through  the  membrane 
produces  an  advection  flux  that  increases  with  current  density.  In 
the  channel,  this  advection  flux  carries  all  species,  including  BH4  as 
shown  in  Fig.  3.  Since  this  model  assumes  only  Na+  and  H2O  pass 
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Fig.  3.  Mole  fluxes  of  BH4  in  the  y-direclion,  half  way  down  the  anode  channel,  as  a 
profile  from  anode  to  membrane.  Positive  fluxes  are  toward  the  membrane  and 
negative  fluxes  are  toward  the  anode.  Results  for  the  baseline  case  at  2.2  V  (peak  power 
density). 

through  the  membrane,  the  net  flux  of  BH4  at  the  membrane 
interface  must  be  zero,  which  can  only  be  satisfied  with  an  elevated 
BH4  concentration  at  the  membrane  interface  to  drive  diffusion 
and  migration  fluxes  which  oppose  the  advection  flux.  This  buildup 
may  also  impact  BH4  crossover,  which  has  been  observed  under 
some  conditions  in  DBFC  experiments  using  Nation  membranes. 
[43], 

Water  crosses  the  membrane  at  a  channel-averaged  rate  of 
8.3  x  10-3  kg  m~2  s-1  in  the  baseline  case  at  2.2  V,  which  is  14.6 
times  larger  than  the  rate  of  water  production  at  the  anode.  Large 
water  crossover  rates  have  been  observed  in  experiments  [44]  and 
can  be  problematic  for  systems  with  large  recirculation  fractions. 
The  permeation  water  flux  (5.5  x  1CT11  kg  itT2  s-1)  is  insignificant 
due  to  the  small  ( ~  5  Pa)  pressure  difference  across  the  membrane. 
The  pressure  difference  across  the  membrane  would  have  to  be 
unpractically  high  for  the  permeation  flux  to  offset  the  electro- 
osmotic  drag  flux. 

Contour  plots  of  species  concentration  in  the  channel  (Fig.  4)  are 
not  a  quantitative  performance  metric,  but  they  can  provide  insight 
into  cell  operation.  For  example,  the  compact  BH4  and  BO2  con¬ 
centration  boundary  layers  near  the  anode  explain  why  the  fuel 
utilization  is  low;  most  fuel  is  carried  through  the  channel  before  it 
can  diffuse  or  migrate  to  the  anode.  Similarly,  the  Na+  concentra¬ 
tion  boundary  layer  at  the  anode  can  be  explained  by  diffusion, 
migration  and  advection  interacting  with  the  electroneutrality 
condition.  The  electric  field  drives  cations  away  from  the  anode,  but 
the  Na+  flux  must  be  zero  there  because  the  wall  is  impermeable 
and  Na+  does  not  participate  in  the  reaction.  The  Na+  concentration 
near  the  anode  is  depressed  to  create  a  diffusion  flux  that  opposes 
the  migration  flux  and  makes  the  net  flux  zero.  Anion  (BH4  and 
BO2 )  concentrations  must  also  decrease  to  maintain  a  neutral  so¬ 
lution,  and  the  lower  anion  concentrations  change  the  rates  of 
transport  and  local  OCV. 

3.2.  Effects  of  reaction  rate  constant 

Fast  forward  reaction  rate  constants  were  chosen  for  the  base¬ 
line  case  to  emphasize  transport  processes,  but  catalysts  in  DBFC 
experiments  are  often  less  active  [15,45],  Slower  anode  rates  were 
examined  to  understand  how  a  less  active  anode  catalyst  would 
affect  cell  performance.  Three  polarization  curves  in  Fig.  5  were 


generated  by  varying  the  anode  forward  reaction  rate  constant  ky, 
between  10°  and  106  m4  mol-1  s_1  with  the  cathode  forward  re¬ 
action  rate  constant  k 2iC  =  10s  m4  mol-1  s_1.  This  range  of  rate 
constants  corresponds  to  a  shift  in  activation  energy  barrier  Eact  of 
approximately  35  kj  mol-1  (0.35  eV),  which  should  encompass  the 
range  of  uncertainty  for  a  selected  catalyst.  The  trend  in  Fig.  5  with 
increasing  reaction  rate  constant  (lower  fact)  indicates  performance 
benefits  that  may  be  obtained  with  the  development  of  more  active 
anode  catalysts  for  the  desired  reactions.  A  faster  anode  reaction 
rate  decreases  the  anode  activation  overpotential,  but  with 
diminishing  returns  for  fcii£I  >10s  m4  mol-1  s_1,  beyond  which  point 
the  anode  activation  overpotential  is  much  smaller  than  other 
overpotentials  in  the  cell. 

A  fourth  polarization  curve  in  Fig.  5  was  generated  using  rate 
constants  from  the  literature  for  BH4  oxidation  on  Pt 
(l<ia  =  10-2  m  s_1  [15])  and  H2O2  reduction  on  Pt 
(/<2,c  =  8  x  10-3  m  s"1  [45]).  This  curve  shows  that  using  Pt  at  both 
electrodes  will  incur  large  activation  losses,  and  peak  power  den¬ 
sity  will  be  lower  (0.175  W  cm-2)  compared  to  the  transport 
dominated  baseline  case  (1.06  W  cm-2).  For  comparison,  power 
densities  of  0.05—0.50  W  cm-2  are  typical  of  DBFC  experiments, 
depending  on  cell  design,  operating  conditions  and  catalyst  selec¬ 
tion  [8], 

Sensitivity  to  reaction  rate  was  evaluated  by  changing  the  for¬ 
ward  reaction  rates  at  both  electrodes  from  the  baseline  values. 
Comparing  ASR  plots  in  Fig.  2(c)  (ky  =  k 2,c  =  106  m4  mol-1  s_1)  to 
Fig.  6  (k^a  =  k 2,c  =  103  m4  mol-1  s~3)  shows  that  the  anode  acti¬ 
vation  ASR  falls  much  faster  than  the  cathode  activation  ASR  when 
the  forward  rate  constants  are  increased.  For  the  conditions 
modeled  here,  the  anode  activation  overpotential  is  therefore  more 
sensitive  to  forward  reaction  rate  constant. 

3.3.  Effects  of  cell  geometry 

For  the  cell  topology  illustrated  in  Fig.  1,  two  dimensions 
pertinent  to  cell  performance  are  channel  and  membrane  thick¬ 
ness.  The  electrolyte  solutions  (in  the  channels)  and  the  membrane 
each  have  an  electrical  resistivity  determined  by  the  mobilities  of 
charged  species  in  the  respective  phase.  The  electrical  resistivity 
incurs  an  ohmic  loss  (and  associated  ASR)  which  is  reduced  as  the 
channel  or  membrane  is  made  narrower.  Furthermore,  a  narrower 
channel  also  results  in  a  steeper  voltage  gradient  in  the  y-direction 
which  produces  a  stronger  electric  field  and  larger  migration  fluxes. 
The  larger  migration  fluxes  can  augment  electrode  reactions  and 
transport  through  the  membrane. 

Fig.  7  plots  power  density  vs.  cell  potential  for  channels  0.50, 
0.75  and  1.00  mm  deep.  The  same  mean  fuel  solution  inlet  velocity 
(0.1  m  s-1)  was  used  in  each  case  to  eliminate  changes  in  the  rate  of 
convection  related  to  velocity.  Shallower  channels  clearly  provide 
higher  power  density  with  little  sacrifice  of  voltage  efficiency.  This 
benefit  is  amplified  further  when  considering  volumetric  power 
density  because  the  shallower  channels  enable  more  cells  to  fit  in  a 
given  volume.  A  further  advantage  of  shallower  channels  is  higher 
channel  fuel  utilization.  With  shallower  channels,  the  concentra¬ 
tion  boundary  layer  occupies  a  larger  fraction  of  the  channel  and 
therefore  less  fuel  flows  through  the  channel  unreacted.  Other 
performance  metrics  are  listed  in  Table  6. 

Most  of  the  ohmic  voltage  loss  between  the  two  electrodes  oc¬ 
curs  in  the  membrane  as  shown  in  Fig.  8,  with  the  losses  being 
highest  near  the  channel  inlet  where  the  current  densities  are 
highest.  This  is  due  to  the  relatively  low  diffusivity  of  Na+  in  Nafion 
(6.93  x  10-10  m2  s-1  vs.  1.33  x  10-9  m2  s-1  in  water  at  infinite 
dilution).  A  thinner  membrane  mitigates  these  losses  and  improves 
the  cell  voltage  and  power  density  as  illustrated  in  Fig.  9,  where 
polarization  and  power  density  curves  are  plotted  for  three 
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Fig.  4.  Contour  plots  of  (a)  BH4,  (b)  BO2  and  (c)  Na+  concentration  (kmol  m  3)  in  the  fuel  channel  for  the  baseline  case  at  2.2  V  (peak  power  density).  Plots  (a)  and  (b)  show  the 
development  of  concentration  boundary  layers  at  the  anode  as  BH4  diffuses  toward  the  anode  and  BO2  diffuses  away  from  the  anode.  In  plot  (c)  the  boundary  layer  near  the  anode 
is  a  result  of  the  zero-flux  boundary  condition  for  Na+,  and  the  boundary  layer  near  the  membrane  is  produced  by  Na+  leaving  the  channel  through  the  membrane. 


membrane  thicknesses.  Reducing  membrane  thickness  from  145  to 
1.45  [tm  (equivalent  to  increasing  membrane  conductivity  by  100 
times)  raises  the  cell  voltage  at  0.40  A  cirT2  from  2.375  V  in  the 
baseline  case  to  2.570  V.  The  very  thin  1.45  um-thick  membrane 
accounts  for  only  3.8— 4.6%  of  the  total  ohmic  losses  in  the  cell 
(depending  distance  from  the  inlets).  Of  course,  the  risk  of  fuel 
crossover  not  considered  in  this  model  would  become  much  more 
significant  with  such  a  thin  membrane  and  would  need  to  be 
included  in  a  final  analysis.  Performance  metrics  as  a  function  of 
membrane  thickness  without  including  the  impact  of  fuel  crossover 
are  listed  in  Table  7. 

3.4.  Effects  of  fuel  concentration  and  flow  rate 

Changes  to  fuel  concentration  and  flow  rate  affect  cell  perfor¬ 
mance  by  changing  concentrations  near  the  electrodes,  thereby 
changing  the  concentration  overpotentials.  Higher  inlet  BH4  con¬ 
centration  provides  a  larger  gradient  to  drive  BH4  from  the  bulk  to 
the  anode,  which  leads  to  higher  fuel  concentration  near  the  anode. 
Higher  flow  rate  improves  convection  transport  of  reactants  to  the 


anode  and  products  from  the  anode.  Both  scenarios  increase  power 
density  by  decreasing  the  concentration  overpotential,  increasing 
current  density,  or  both. 

A  plot  of  power  density  vs.  cell  voltage  (see  Fig.  10)  for  a  range  of 
reasonable  BH4  inlet  concentrations  between  0.1  and  0.5  M  shows 
that  peak  power  density  is  approximately  proportional  to  the  BHT 
inlet  concentration.  The  rate  of  peak  power  density  increase  with 
concentration  slows  at  higher  concentrations  (see  Fig.  10  inset) 
because  other  losses  in  the  cell  (ohmic  and  concentration  over¬ 
potentials)  play  a  larger  role  with  increasing  current  density. 
Because  the  inlet  BH4  concentration  affects  cell  performance  by 
altering  the  concentration  overpotential,  the  strength  of  the  effect 
depends  on  cell  current  density.  At  high  current  density  the  con¬ 
centration  overpotential  is  large,  and  increasing  the  BHf  inlet 
concentration  can  increase  the  power  density  substantially.  The 
fuel  utilization  at  high  current  density  does  not  change  with  inlet 
BH4  concentration  because  both  the  rate  of  BH4  flow  through  the 
cell  and  the  rate  of  consumption  are  proportional  to  concentration. 
At  low  current  density  the  concentration  overpotential  is  already 
small  compared  to  other  losses  in  the  cell,  so  there  is  little 
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Fig.  5.  Polarization  curves  with  anode  forward  reaction  rate  constants  of  10°,  103  and 
106  m4  mol-1  s  1  and  a  cathode  forward  reaction  rate  constant  of  106  m4  mor1  s 
The  bottom  curve  assumes  forward  reaction  rate  constants  from  Finkelstein  et  al.  on 
Pt;  anode:  10  2  m  s  1  and  cathode:  8  x  10  3  m  s  All  parameters,  other  than  rate 
constants,  are  the  same  as  in  the  baseline  case,  given  in  Table  4. 


Cell  Voltage  [  V  ] 

Fig.  7.  Average  power  density  vs.  cell  potential  for  channels  that  are  0.50,  0.75  and 
1.00  mm  deep.  Aside  from  channel  depth,  all  parameters  are  the  same  as  in  the 
baseline  case,  which  is  labeled  "BL".  The  bold  dashed  curve  is  drawn  through  the 
maximum  power  density  at  each  channel  depth  and  describes  the  combinations  of 
channel  depth  and  cell  voltage  which  yield  peak  power.  In  each  case  the  mean  velocity 
of  fuel  solution  entering  the  fuel  channel  was  0.1  m  s 


opportunity  for  higher  BH4  inlet  concentration  to  affect  perfor¬ 
mance.  Furthermore,  at  low  current  density  the  additional  fuel  goes 
unreacted,  substantially  decreasing  fuel  utilization. 

The  choice  of  inlet  fuel  concentration  in  a  DBFC  system  depends 
on  whether  a  recirculation  loop  is  used.  Recirculation  is  necessary 
in  practical  DBFCs  to  achieve  acceptable  system-level  fuel  utiliza¬ 
tion.  Given  practical  lower  bounds  on  fuel  flow  rate  and  channel 
depth,  the  single-pass  fuel  utilization  will  be  low  —  in  the  case  of 
this  study,  no  greater  than  10%.  If  a  recirculation  loop  is  used,  it  may 
be  preferable  to  store  the  fuel  at  high  concentration  near  the 
saturation  limit  and  add  it  to  the  recirculation  loop  slowly  to 
maintain  low  (~0.1  M)  concentration  in  the  cell.  The  low  fuel 
concentration  in  the  cell  will  yield  lower  peak  power  density  (as 
shown  in  Fig.  10)  but  will  mitigate  losses  to  hydrolysis  (Reaction 

(R4)). 

Power  density  is  plotted  vs.  cell  voltage  for  three  fuel  flow  rates 
in  Fig.  11.  The  power  density  increases  with  increasing  flow  rate  at 

10’  -1 - 


io°  - 


-  Anode  Activation  -  Cathode  Activation 


Anode  Concentration - Cathode  Concentration 

^q.3  _ - Fuel  Ohmic  - Oxidizer  Ohmic 

- Membrane  Ohmic 

0  50  100  150  200 

Distance  from  Inlet  [  mm  ] 

Fig.  6.  Area  specific  resistance  at  points  down  the  channel,  for  the  case  with  forward 
reaction  rates  set  to  103  and  cell  voltage  of  2.2  V.  Aside  from  the  reaction  rate  con¬ 
stants,  all  parameters  are  the  same  as  in  the  baseline  case. 


all  current  densities  as  the  higher  rates  of  transport  change  con¬ 
centrations  near  the  anode.  The  pressure  drop  in  the  channel  in¬ 
creases  in  proportion  to  the  flow  rate,  but  the  fuel  utilization  falls  as 
indicated  in  Table  8.  Increasing  the  fuel  flow  rate  by  a  factor  of  four 
increases  the  pressure  drop  (pumping  losses)  by  a  factor  of  four,  but 
only  increases  the  peak  power  by  48%. 

3.5.  Implications  for  cell  design 

The  results  in  Sections  3.1— 3.4  can  guide  DBFC  design  by 
showing  how  performance  varies  with  operating  conditions  and 
cell  design  parameters.  For  conditions  in  which  cell  performance  is 
controlled  by  transport,  the  assumptions  regarding  reaction  ki¬ 
netics  and  OCV  are  reasonable,  and  the  model  can  be  used  to 
compare  specific  design  options.  For  example,  Fig.  10  can  be  used  to 
choose  an  inlet  BH4  concentration  by  drawing  a  line  through  the 
peak  power  density  of  each  curve.  The  dashed  line  in  Fig.  10  de¬ 
marcates  the  trade  space  maximizing  power  density  as  a  function  of 
BH4  inlet  concentration  and  cell  voltage.  The  region  to  the  right  of 
the  dashed  line  describes  designs  or  operating  points  which  trade 
power  density  for  higher  voltage  efficiency.  Since  the  energy  den¬ 
sity  of  a  system  is  proportional  to  fuel  concentration,  voltage  effi¬ 
ciency  and  fuel  utilization,  moving  to  the  right  in  Fig.  10  can 
substantially  increase  the  system  energy  density.  The  maximum 
power  density  achievable  given  an  energy  density  requirement  (or 
vice  versa)  can  be  estimated  from  Fig.  10. 

In  general,  shallower  channels  are  preferable  because  they  in¬ 
crease  system-level  power  and  energy  density.  Power  density  is 
improved  when  shallower  channels  decrease  the  stack  volume,  and 


Table  6 

Performance  metrics  with  respect  to  channel  depth.  All  other  parameters  are  same 
as  the  baseline  case. 


Parameter 

Channel  depth  [mm] 

0.50 

0.75 

1.00 

Peak  power  density  [W  cm-2] 

1.28 

1.15 

1.06 

Cell  voltage@peak  power  density  [V] 

2.16 

2.19 

2.21 

Fuel  utilization@peak  power  density  [%] 

10.2 

6.05 

4.14 

Power  density@2.5  V  [W  cm-2] 

0.82 

0.77 

0.73 
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Fig.  8.  Electric  potential  profiles  across  the  cell  from  anode  to  cathode,  near  the  inlet 
and  near  the  outlet.  The  membrane  clearly  makes  the  largest  contribution  to  the  total 
ohmic  loses,  and  the  decrease  in  ohmic  losses  from  inlet  to  outlet  is  related  to  the 
decrease  in  current  density.  Data  in  this  figure  are  from  the  baseline  at  2.2  V. 

energy  density  is  improved  when  the  smaller  stack  frees  additional 
volume  for  reactant  storage.  When  the  additional  volume  is  used  to 
add  active  area  to  the  stack,  then  it  can  produce  the  same  power  at 
lower  current  density  and  operate  more  efficiently  (higher  voltage 
per  cell).  These  trends  are  seen  in  the  DBFC  model  results.  For 
example,  a  stack  with  channels  1.00  mm  deep  operating  at  2.21  V 
per  cell  (point  A  in  Fig.  7)  would  have  the  same  volumetric  power 
density  as  a  stack  with  twice  as  many  0.50  mm  deep  channels 
operating  at  2.64  V  per  cell  (point  B  in  Fig.  7),  but  the  stack  with 
shallower  channels  would  have  36%  less  voltage  loss.  Channel  fuel 
utilization  is  lower  in  the  shallow  channel  case  (3.2%  vs.  4.2%) 
because  of  the  lower  current  density,  so  a  complete  system  trade 
would  have  to  include  the  recirculation  fraction  and  fuel  storage 
density  to  find  an  optimum  channel  depth. 

Another  performance  trade  is  the  increase  in  peak  power  vs. 
parasitic  power  consumed  by  the  fuel  and  oxidizer  pumps.  Table  8 
shows  that  pressure  drop  through  the  fuel  channel  is  directly 


Table  7 

Performance  metrics  with  respect  to  membrane  thickness.  All  other  parameters  are 
same  as  the  baseline  case. 


Parameter 

Membrane  thickness  [pm] 

1.45 

35.25 

72.50 

108.75 

145.00 

Peak  power 

1.21 

1.17 

1.13 

1.10 

1.06 

density  [W  cm-2] 

Cell  voltage@peak  power 

2.43 

2.37 

2.32 

2.26 

2.21 

density  [V] 

Fuel  utilization@peak 

4.26 

4.24 

4.19 

4.17 

4.14 

power  density  [%] 

Power  density@2.5  V  [W  cm-2] 

1.18 

1.07 

0.94 

0.82 

0.73 

proportional  to  the  fuel  flow  rate,  as  expected  for  laminar  incom¬ 
pressible  flow  through  a  straight  channel.  The  peak  power  is  also 
proportional  to  the  fuel  flow  rate,  but  increases  at  a  rate  of  less  than 
8  mW  crrr2  per  mL  miiT1  increase  in  fuel  flow  rate.  Thus,  there 
may  be  a  design  point  beyond  which  the  additional  power  from 
faster  flow  is  not  worth  the  increase  in  pressure  drop  and  loss  of 
fuel  utilization. 

The  low  channel  fuel  utilization  in  all  cases  in  this  study  bolsters 
experimental  observations  that  any  practical  DBFC  will  require 
recirculation  of  the  reactants  to  achieve  high  (>90%)  overall  fuel 
utilization.  The  near  irreversibility  of  the  borohydride  oxidation 
reaction  makes  the  reaction  rate  insensitive  to  BOJ  concentration; 
reaction  rates  varied  by  less  than  0.3%  in  simulations  with  BO2  inlet 
concentrations  between  1CT6  M  and  0.3  M.  The  insensitivity  to  BO2 
concentration  makes  extrapolating  from  once-through  simulation 
results  to  a  recirculating  system  straightforward  because  the  only 
parameter  related  to  both  cell  performance  and  the  recirculation 
fraction  is  the  BH4  concentration  at  the  inlet.  A  simple  recirculation 
model  was  constructed  which  interpolates  among  the  results  from 
the  DBFC  model  with  varying  inlet  concentration  of  BH4.  The 
recirculation  model  was  used  to  investigate  the  relationships  be¬ 
tween  anode  recirculation  volume  fraction,  fuel  utilization  and 
power  density.  The  results  for  one  case  are  shown  in  Fig.  12,  in 
which  the  NaBFLi  fuel  is  stored  at  0.5  M,  the  cell  voltage  is  2.5  V  and 
fuel  solution  is  added  upstream  of  the  channel  at  the  same  volu¬ 
metric  rate  waste  is  rejected  downstream  of  the  channel.  Fig.  12 


Fig.  9.  Polarization  and  power  density  curves  for  membrane  thicknesses  of  1.45,  72.5, 
and  145  pm.  All  other  parameters  are  the  same  as  in  the  baseline  case,  which  is  labeled 
“BL”.  The  thickness  of  fully  hydrated  Nafion  115  in  the  Na+  form  is  145  pm  [52],  Note 
that  decreasing  the  membrane  thickness  is  approximately  equivalent  to  increasing  the 
membrane  conductivity  by  the  same  factor. 


Fig.  10.  Average  power  density  vs.  cell  potential  for  fuel  inlet  concentrations  ranging 
from  0.1  M  to  0.5  M.  All  other  parameters  are  the  same  as  in  the  baseline  case,  which  is 
labeled  “BL".  The  bold  dashed  curve  is  drawn  through  the  maximum  power  density  at 
each  concentration  and  describes  the  combinations  of  inlet  concentration  and  cell 
voltage  which  yield  peak  power. 
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Fig.  11.  Power  density  vs.  cell  voltage  curves  for  fuel  solution  flow  rates  of  15,  30  and 
60  mL  min-1.  All  other  parameters  are  the  same  as  in  the  baseline  case,  which  is 
labeled  “BL''.  The  bold  dashed  curve  is  drawn  through  the  maximum  power  density  at 
each  flow  rate  and  describes  the  combinations  of  flow  rate  and  cell  voltage  which  yield 
peak  power. 


Table  8 

Performance  metrics  with  respect  to  fuel  flow  rate.  All  other  parameters  are  same  as 
the  baseline  case. 


Parameter 

Inlet  fuel  flow  rate  [mL  min  *] 

15 

30 

60 

Peak  power  density  [W  cm-2] 

0.72 

0.88 

1.06 

Cell  voltage@peak  power  density  [V] 

2.37 

2.29 

2.21 

Fuel  utilization@peak  power  density  [%] 

10.5 

6.60 

4.14 

Power  density@2.5  V  [W  cm-2] 

0.65 

0.70 

0.73 

Channel  pressure  drop  [Pa] 

51.7 

104 

210 

shows  that  the  fuel  utilization  improves  by  a  factor  of  more  than  5 
over  a  once-through  approach  when  99%  of  the  fuel  flow  is  recir¬ 
culated,  at  the  cost  of  a  50%  decrease  in  power  density  due  to  the 
lower  inlet  BH4  concentration. 

4.  Conclusions 

A  2-D  DBFC  model  has  been  developed  which  captures  the 
major  losses  (activation,  ohmic  and  concentration)  and  down-the- 
channel  effects  in  a  DBFC  with  simple  cell  geometry.  The  results 
presented  here  assume  simple  electrode  reaction  models,  but  the 
model  framework  can  accommodate  competing  reactions  and 
mixed  potentials  when  the  mechanism  and  rate  parameters  for 
borohydride  oxidation  are  better  understood.  The  DBFC  model  was 
used  to  explore  performance  trends  for  conditions  in  which  the 
ideal  reactions  are  good  approximations  to  the  actual  rates.  The 
results  show  that  transport  due  to  advection,  diffusion  and 
migration  all  play  an  important  role  in  the  boundary  layers  devel¬ 
oped  at  the  electrode-channel  and  membrane— channel  interfaces. 
The  relative  importance  of  these  transport  mechanisms  evolves  as 
the  flow  proceeds  down-the-channel  because  concentration 
boundary  layers  develop  and  strongly  affect  cell  behavior  by 
shifting  local  concentrations  and  equilibrium.  Models  attempting  to 
predict  DBFC  performance  in  cells  with  geometry  similar  to  Fig.  1 
must  account  for  the  change  in  conditions  at  points  far  from  the 
inlet,  preferably  in  an  explicit  manner  as  in  this  work. 

Cells  similar  to  Fig.  1  should  have  channels  as  shallow  as  prac¬ 
ticable  to  increase  fuel  utilization  and  decrease  ohmic  losses  in  the 


Recirculation  Fraction 

Fig.  12.  Channel  fuel  utilization,  overall  fuel  utilization,  BHt  and  BO2  concentrations  at 
the  inlet,  all  as  functions  of  the  fuel  recirculation  volume  fraction.  The  concentration  of 
BH4  in  the  fuel  added  to  the  recirculation  loop  was  assumed  to  be  0.5  M.  Aside  from 
the  inlet  concentrations,  all  parameters  are  the  same  as  in  the  baseline  case. 

fuel  and  oxidizer  solutions.  DBFC  performance  (peak  power  density, 
voltage  efficiency  and  fuel  utilization)  rises  more  quickly  than 
pumping  losses,  so  any  channel  depth  reduction  yields  a  net  gain  in 
performance.  Membranes  are  similar  in  that  thinner  membranes 
are  preferable  in  the  absence  of  fuel  crossover.  Higher  fuel  flow 
rates  yield  higher  peak  power  density  with  slightly  improved 
voltage  efficiency,  but  with  lower  fuel  utilization  and  increased 
pumping  losses  which  vary  with  cell  geometry.  Water  crossover 
from  fuel  channel  to  oxidizer  channel  must  be  addressed  in  a 
practical  DBFC  system  using  recirculation. 

The  sensitivity  of  DBFC  performance  to  design  and  operating 
parameters  encourages  the  use  of  a  model  when  designing  for  a 
specific  application.  As  illustrated  here,  there  are  many  ways  to 
adjust  the  power  and  similarly,  energy  density  of  a  DBFC  system  by 
trading  peak  power  density,  voltage  efficiency,  fuel  utilization,  and 
recirculation  fraction.  A  model  can  guide  DBFC  designers  to  a  point 
closer  to  the  final  design  before  starting  expensive  prototype 
design-build-test  iterations.  Future  work  to  calibrate  the  model 
against  DBFC  experiments  and  to  incorporate  a  more  comprehen¬ 
sive  multi-step  reaction  mechanism  into  the  model  will  provide 
increased  confidence  that  the  calibrated  model  can  be  used  to 
design  DBFC  systems  for  various  applications. 
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Appendix  A.  Nomenclature 

a  anodic  parameter  (subscript) 

/?  electrochemical  reaction  symmetry  factor 

c  cathodic  parameter  (subscript) 

ch  channel  discretization  (subscript) 

Cfc  concentration  of  species  k 

5k  apparent  molar  volume  of  species  k  in  water 

D/(  diffusivity  of  species  k 

standard  half  cell  potential  for  reaction  R 
e_i  electrode  interface  discretization  (subscript) 
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F  Faraday’s  constant 

/  FIRigT 

gk  reaction  rate  order  with  respect  to  species  k 

h  channel  depth 

i  current  density 

J  mass  flux 

J  k  mole  flux  of  species  k 

kq  rate  constant  for  reaction  q 

5  roughness  factor 

p  kinematic  viscosity 

m_i  membrane  interface  discretization  (subscript) 

rid  electro-osmotic  drag  coefficient 

ne  number  of  electrons  in  the  rate  limiting  reaction  step 

ne iec  number  of  electrons  transferred  per  mole  of  fuel 

consumed 

<p  electric  potential 

P  pressure 

Ph2o  membrane  permeability  to  water 

p  solution  mass  density 

pc  net  charge  density 

R  residual  (error)  in  model  governing  equation 

Rig  ideal  gas  constant 

rq  rate  of  reaction  q 

T  temperature 

uk  mobility  of  species  k 

v  velocity 

Vfc  stoichiometric  coefficient  of  species  k 

V  volumetric  flow  rate 

Wk  molar  mass  of  species  k 

x  coordinate  axis  parallel  to  membrane 

y  coordinate  axis  perpendicular  to  membrane 

Yk  mass  fraction  of  species  k 

zk  charge  number  of  species  k 
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